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We discuss an ensemble of topological solitons - instanton-dyons and antidyons - in SU{2) pure 
gauge theory at finite temperatures above and below the deconfinement phase transition tempera- 
ture. The main focus is on the combined effect of this ensemble on the so called effective holonomy 
potential, which drives the confinement/deconfinement phase transition. Using a simple model with 
excluded volume and lattice data on caloron density we find that repulsive part of the potential is 
robust enough to induce the phase transition at the right temperature. Model's predictions - the 
holonomy potential, electric and magnetic screening masses as a function of T - are in qualitative 
agreement with the available lattice data. Further predictions are densities of various dyon types as 
a function of temperature: while some lattice measurements of them had been made, much more 
accurate data are needed to test these predictions. 

PACS numbers: 



INTRODUCTION 

The expectation value of the Polyakov loop P = 
Pexp[z ^ Aqg^t] is an order parameter of confine- 
ment/deconfinement phase transition in pure Yang Mill 
theory. The gauge invariant integral over the thermal 
(Matsubara) circle is also called holonomy. Its average 
value and fluctuations are studied numerically in lattice 
simulations, or described by an effective potential U{P) 
fitted to lattice data, see e.g. [ ]. This potential plays 
a prominent role in current models of finite tempera- 
ture QCD, such as Polyakov-Nambu-Jona-Lasinio model 
[2, 3], yet the understanding of its origin is still miss- 
ing. Some recent attempts to understand the structure 
of this potential have been made, both using the lattice 
gauge theory and using functional renormalisation group 
method [ - ]. 

Instantons, the 4-dimensional solitons carrying the 
topological charge, are known to be important compo- 
nent of the gauge field theory. They have been used suc- 
cessfully in the description of many nonperturbative phe- 
nomena (for review see ["]), including U{1) and SU{Nf) 
chiral symmetry breaking, but not confinement. At 
nonzero asymptotic holonomy Ao{r oo) 7^ the in- 
stanton solution is modified and as a result the instanton 
is split into Nc constituents known as KvBLL self-dual 
instanton-dyons [ , ]. 

Interactions and statistical ensemble of instanton- 
dyons and antidyons have been studied qualitatively us- 
ing neutral dyon-antidyon pairs in [ ], recently comple- 
mented by first numerical simulations [ ]. These works 
used holonomy as an input, defining masses and other 
properties of the instanton-dyons. They had studied the 
interaction between self-dual and anti-self- dual sectors 
via light fermions, and focused on the fermion collec- 
tivization and chiral symmetry breaking phase transition. 
In this paper we turn our attention to the back reaction 



of the topology on the holonomy potential, via a simple 
model with excluded volume, much in a spirit of van der 
Waals theory of non-ideal gases. 

Two distinct sectors of instanton interaction are 
treated quite differently. Pure sefl-dual (or anti-self- 
dual) sector has classicaly degenerate moduli spaces pa- 
rameterized by collective coordinates: their studies are 
rather complete, and their nontrivial geometry were ex- 
tensively studied in mathematical literature since the 
1970s by Atiyah, Hitchin and others. 

The instanton-anti-instanton (//) configurations are 
studied significantly less. The moduli spaces are in this 
case substituted by the "streamline configurations" [12, 

] which smoothly interpolate between the separated 
instanton-anti-instanton pair and the perturbative fields. 
Close // pairs correspond to weak fields, which cannot 
be treated semiclassically and should be subtracted from 
the semiclassical configurations. This physical idea has 
been implemented in the Instanton Liquid Model via an 
"excluded volume" , which generates a repulsive core and 
stabilizes the density. 

In a few important cases, in which the partition func- 
tion is independently known, such subtraction can be 
performed exactly, without any parameters. The // pair 
contribution to the partition function in QM instanton 
problem has been done via the analytic continuation in 
the coupling constant —g^ by Bogomolny [i-t] and 

Zinn- Justin [ ] (BZJ), who verified it via known semi- 
classical series. Another analytic continuation has been 
used by Balitsky and Yung [ ] for supersymmetric quan- 
tum mechanics. 

Recently Poppitz, Schafer and Unsal (PSU) [16, 17] 
used BZJ approach in the N = 1 Super- Yang-Mills 
theory on x S^^ observing that the result obtained 
matches exactly the result derived via supersymmetry 
[ ]. PSU papers are the most relevant for this work, 
as they focus on the instanton-dyons (referred to as 
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instanton-monopoles in their work). The issue in this 
case is more comphcated than for //: dyon-antidyon pair 
always has uncompensated charges, magnetic or electric. 
Studies of the corresponding "streamline configurations" 
is yet to be done. 

In a carefully tuned weakly coupled, supersymmetric 
setting, softly broken by a small gluino mass, PSU [i6, 17] 
have calculated the contribution of instanton-dyons and 
neutral instanton-dyon pairs. (They call them neutral 
bions^ but they are also known as LL, MM molecules in 
other works.) They observed that the repulsion in the 
latter term can overcome the former one and drive the 
confinement /deconfinement phase transition. 

Supersymmetry forces the perturbative holonomy po- 
tential to be canceled, between contributions of gluons 
and (periodic) gluinos. Pure Yang-Mills (we are inter- 
ested in) has no gluinos and there is no cancellations of 
the perturbative potential, which prefers the deconfined 
phase. The QCD-like theories with thermally compact- 
ified fermion (i.e. fermions with anti-periodic temporal 
boundary condition) further add to this perturbative de- 
confining potential. In order to generate confinement, 
one needs to find some nonperturbative mechanism which 
is strong enough to compete with the perturbative holon- 
omy potential. 

In this work we show that the free energy of the 
instanton-dyons can induce holonomy potential which 
produces confinement in the pure gauge theory (us- 
ing as an example the simplest SU(2) case) in qualita- 
tive agreement with the lattice data. The central role 
is played by the effective repulsive interaction between 
dyon-antidyon, as in the PSU works. We should stress, 
however, that PSU in [ ] did consider non supersymet- 
ric, pure Yang-Mills case. We comment on the similar- 
ities and the differences of their work and ours later in 
the text. 

For definiteness, we start with fixing global parameters 
of the SU (2) gauge theory topology, using lattice works 
[19] and [20]. The finite-T instantons or calorons possess 
only the topological charge - they have neither electric 
nor magnetic charges and thus no direct coupling to ei- 
ther electric or magnetic holonomies 6, a (see next sec- 
tion) . They can be identified on the lattice via a number 
of well developed methods, e.g. the so called "cooling" . 
In Fig.l we present the data from [ ] and their fit by the 
semiclassical expression for the dimensionless density (in 
units of T^) of instantons/calorons 

22 /T\ 

with parameters^ K = 0.024, A/Tc = .36. The caloron 
action at Tc is 7.50, so per dyon it makes Sd = Scaif^ — 



Lattice practitioners usually fix "physical units" via T = string 




FIG. 1: Caloron density as a function of T/Tc. The solid 
curve is the semiclassical fit ricai = KScai^'^'"'^^ in units of T 
with parameters K = 0.024, Seal = Stt'^ / g'^ (T) , open (filled) 
points are the lattice data from [19] ([20]). 



3.75, which gives an idea how semiclassical the discussed 
objects are. (SU(3) instantons have actions Seal ^ 12 or 
Sd = Scai/^ ^ 4, quite close in magnitude.) After those 
parameters are fixed, one knows semiclassical densities of 
the dyons and their pairs, as we explain in detail below. 



THE HOLONOMY POTENTIAL 

In two-color SU{2) gauge theory there is only one 
diagonal color matrix, thus one electric and one mag- 
netic holonomy. An electric holonomy is closely related 
to the local Polyakov loop, which for SU(2) group can 
be parametrized by a single (space dependent) angle 

e = v/T 

2TrP(.)=cos(J^j , ^=^(^-1) (2) 

where P{x) = V exp [i f AqcIt) . The re-definition of v 
into field 6, by a shift and re-scaling, is done for fur- 
ther convenience (note that our 6-field is labeled b' in 
[16]). Zero v means trivial (deconfining) holonomy, zero 
b means confinement. 

Magnetic potential (or magnetic holonomy) a was in- 
troduced by Polyakov Ref.[ ], together with an obser- 



tension, to which the QCD value is ascribed. For pure gauge 
SU(2) one then finds Tc ^ 270 MeV. 
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= J Da exp —{dcr)'^/2-\-ieicr{x)-\-ie2C 
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which expresses Coulomb forces via corresponding po- 
tential field (J. Since our objects have both electric and 
magnetic charges and potentials, the trick is used twice. 
The amplitudes for the 4 types of dyons take the expo- 
nent al form [16] 



J[4 ^-b+ia-Sd 
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Note that there is no i in the electric part. This be- 
havior is not however because of the Euclidean formu- 
lation of the theory, as the real charge couples the elec- 
tric potential Aq as exp(i J A^dr) both in Euclidean and 
Minkowski space-time (see the discussion in ["^']). In fact 
the coupling to 6-field is through the action of instanton- 
dyons which is h dependent, and it always appears in the 
combination Sd which is the action of the L and M 
dyon. Contributions of all four instanton-dyons can be 
combined into the following effective potential 



eff 



-4n2 cos(cr) cosh(6) 



(6) 



where we use a "generic instanton-dyon density" = 
i'Ticai)^^^'' as a prefactor. The Coulomb gas partition 
function is recovered if 6, a are integrated out with a ki- 
netic energy Ck = 2(4^ [i^icr)'^ + (9^6)^] . 

The next step is an introduction of the binary repul- 
sion and the excluded volume in LL and MM channels. 
Subtracting pair contraction of instanton-dyon-antidyon 
the partition function takes the form 
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where A is the subtracted amplitude 



A = 47r 



dr r^e 



-Vddir) 



(7) 



(8) 



where the upper sign refers to M type and the lower sign 
to L type instanton-dyons and anti-dyons and. The effec- 
tive dyon-antidyon interaction is Coulombic at large dis- 
tances Vdd = —Sirp/g'^r and regulated at small distances 
by some particular scheme, e.g. 1/r + 
needed to stabilise the coulomb plasma. If the objects 
were for example atoms this regularization would suffice. 



It is the instantonic nature of these objects - not regu- 
larisation - which requires the additional subtraction as 
the close pair of instanton-dyon and antidyon are pertur- 
bative fields, and it is this subtraction which brings in a 
new parameter tq or A. 

The effective potential then becomes 

Veff = 2nd{—2cosacoshb-\-ndAcosh.{2b)) . (9) 



and has a minimum at 



0, brniyi 




for 2ndA > 1 
for 2ndA < 1 



(10) 

This is the key point: confining regime (the one with 
bmin = 0) can only be reached if the instanton-dyon den- 
sity is sufficiently large. 

The third step is inclusion of the one-loop perturbative 
potential [ .z] 



pert 
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(11) 



which by itself prefers a trivial holonomy v = O^b = Sd- 
Including this terms one finds that the phase transition is 
pushed to higher densities. By demanding that the sec- 
ond order coefficient of b vanishes, we obtain the equation 
for critical temperature of confinement/deconfinement 
transition 



c\2 



= 0, 



(12) 



Note that all (dimensionfull) quantities above are in 
units of T. Since we take that = \/ncai(Tc) = 
^fKj2Scai{Tcfe-^^-'^^-^ ^ 0.145 we need to have A ^ 
5. While A is much larger than the dyon volume ~ 
47r/3(l/7r)^, including the Coulomb enhancement factor 
~ exp(0(6'c^)) expected for the pair of "half- annihilated" 
dyon-antidyon gives this parameter the correct order of 
magnitude. 

Let us improve the model a bit by noting that the ex- 
cluded volume should be different for M and L dyons, 
and in fact should depend on the holonomy 6, in such 
a way that ro cx {1 ± b/Sd)~^' To account for this 
^(1 ± where ± refers to the instanton-dyon with 

factors e^^ respectively. Further each dyon carries with it 
a moduli space metric parameter in front of the densities, 
which should also be included and that can be obtained 
by replacing ride^^ rid{l ± ^)e^^. This results in a 
somewhat longer expression for the effective potential 
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(13) 
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The effect of these factors helps the phase transition to 
occur and similar analysis as before gives A ~ 2.3. The 
resulting holonomy potential at T/Tc = 0.8,1., 1.5 are 
shown in the upper plot of Fig. 2. 

To conclude this section, let us mention that the po- 
tentials (9) and (13) are similar to the one discussed in 
[16] for the case of pure Yang-Mills. We briefly comment 
on two important differences: the "excluded volume" pa- 
rameter A and the absence of the magnetic bion cos(2cr). 
The magnetic bion term in our analysis can be obtained 
in a similar way to the subtraction, by a virial-type of 
expansion. However since the ML and LM pair repel, 
this term will be accompanied by a coulomb suppres- 
sion term, rather then coulomb enhancement. Therefore 
we believe that such terms in the effective action can 
be safely neglected. We also note that the coefficient 
A should formally be temperature dependent, at least 
through the coupling constant. The naive way to incor- 
porate this dependence would seem to suggest that this 
term is growing with temperature, which is unphysical, 
as it is expected that instanton-dyon effects and charges 
become screened and the effects on the holonomy be- 
come less and less important. Since our paper focuses 
at describing the physics around Tc, we assumed the pa- 
rameter A to be constant. 



PREDICTIONS 

Now that the model is fully defined, some of its predic- 
tions can be compared to the lattice data. We start with 
two important observables usually not associated with 
topology and measured directly from propagators: the 
(mean-field, bare) electric and magnetic screening masses 



eff 



962 



b^brr 



(14) 
(15) 



b=brr 



where bmin^ ^min are values where the potential takes its 
minimum value. 

Their comparison with lattice data is shown in the 
lower plot of Fig. 2. Our model produces phase transition 
of the second order, as is known to be the case for SU(2) 
theory. This implies that the bare electric mass must 
vanish at Tc. While this is indeed evident from holonomy 
distributions observed on the lattice (see e.g. Fig. 4 of 
[19]), the electric mass obtained from gluon propagators 
[26] (shown by the open squares in Fig. 2) indicate only 
partial downward shift. Note, however, that the lattice 
measurements are for the full, physical mass of the glu- 
onic propagator, and that our model has infinitely many 
6-field self-couplings which should renormalize the mass 
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FIG. 2: The upper plot shows the effective potential Veft{h) /T 
(13) for T/Tc = 0.8, 1, 1.5 shown by the dashed,solid and dot- 
dashed lines, respectively. The plot shows electric ttie/T and 
magnetic vtlmIT screening masses versus temperature, indi- 
cated by the solid and dashed lines, respectively. Thick lines 
are our model, the data points are from lattice propagators 
[26], the lines connecting data points are shown simply for 
their identification. 



especially around where the fluctuations are most im- 
portant (and create critical indices different from mean 
field ones, as e.g. observed on the lattice in [ ]). 

The behavior of the magnetic mass in the model has 
a fairly smooth behavior through the phase transition, 
with a small kink at Tc which is due to the dependence 
of the magnetic mass on the density the M and L dyons 
which exhibits a non smooth behavior due to the 2nd or- 
der phase transition. It can be seen that the SU(2) data 
on magnetic mass [ ], shown by triangles, have only 
two points at two temperatures. However analogous but 
much more complete SU(3) data shows a smooth mM(T) 
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FIG. 3: Prediction of the model for the temperature depen- 
dence of the density of the instanton-dyons are shown by the 
hnes, those with sohd and dashed Unes are for M, L type 
dyons, respectively. Open (filled) circles show identified M- 
type dyons from ref. [ 9] ([20]). The crosses show "unidenti- 
fied topological objects" from [19]. Circles and crosses provide 
the lower and the upper bound for the dyon density. 



behavior, in spite of the 1-sf order transition in this the- 
ory. Note that the absolute scale of the masses are about 
factor 2 off, which is hardly surprising since the normal- 
ization was done in a topological sector with similar un- 
certainty (see Fig.l). The ratio ttim/'^e predicted by 
the model and from the lattice agree rather well. We 
thus conclude that the model performs at a qualitative 
level well enough. 

More specific predictions of the model are of course the 
dimensionless (i.e. in units of T) densities of monopole- 
dyons of particular types: 

riM^L = nd{l ± bmin/Sd)e^^^'- 

X (^1 - ^Audil ± brr^in/Sd)-^e^''--^ (16) 

The densities are shown as curves in Fig. 3. At T > Tc 
there are two curves, identifying lighter (and thus more 
numerous) M dyons as well as L ones. Below Tc the 
curves collapse into one single curve as the potential de- 
velops a minimum at 6 = 0. 

Two examples of the best lattice efforts to determine 
the densities in SU(2) pure Yang-Mills theory are refs. 
[19] and [20]. They use 24^ x 6 and 20^ x 4 lattices, the 
former with couplings tuned to three points with T < Tc 
and two T > Tc. The L and M dyons can in principle 
be identified by their electric and magnetic charges as 
well as by the value of the Polyakov loop at their centers 
taking values ±1 for the two type of dyons. In practice. 



the heavy L dyon has such a small size that it cannot 
be identified by this method at current lattices. The 
M instanton-dyons, however, were identified and the re- 
sults are shown as circles in Fig. 3. The authors of [19] 
note however that the efficiency of M-dyon identification 
is quite low, and depends on the temperature. Crosses 
show "unidentified topological objects" from the same 
work [ ], to be taken as an upper bound. As one can 
see from the plot, our predictions lie between the lower 
and upper bound. More work is obviously needed to 
test these predictions. Apart from the densities of the 
instanton-dyons, quite valuable would be lattice studies 
of their spatial correlations, as those can help understand 
their interactions as well. 

Concluding this work, let us comment that while the 
model proposed is rather schematic, it represents a poten- 
tially important link between the gauge topological sector 
- instanton-dyons - and confinement in QCD-like theo- 
ries. Its more quantitative forms and generalization to 
QCD with fermions of different kind seem to be straight- 
forward. Adding fundamental fermions would induce ex- 
tra binding of LL pairs via fermonic zero modes: this 
will shift the deconfinement transition to lower T and/or 
higher instanton-dyon densities, as is indeed observed 
phenomenologically. We intend to study topology and 
all related phenomena quantitatively, by a direct Monte- 
Carlo simulations, elsewhere. 
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